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We investigate QCD-like theory with exact center symmetry, with emphasis on the 
finite-temperature phase transition concerning center and chiral symmetries. On the 
lattice, we formulate center symmetric SU{3) gauge theory with three fundamental 
Wilson quarks by twisting quark boundary conditions in a compact direction (Z 3 - 
QCD model). We calculate the expectation value of Polyakov loop and the chiral 
condensate as a function of temperature on 16^ x 4 and 20^ x 4 lattices along the 
line of constant physics realizing mps/my = 0.70. We find out the first-order cen¬ 
ter phase transition, where the hysteresis of the magnitude of Polyakov loop exists 
depending on thermalization processes. We show that chiral condensate decreases 
around the critical temperature in a similar way to that of the standard three-flavor 
QCD, as it has the hysteresis in the same range as that of Polyakov loop. We also 
show that the flavor symmetry breaking due to the twisted boundary condition gets 
qualitatively manifest in the high-temperature phase. These results are consistent 
with the predictions based on the chiral effective model in the literature. Our ap¬ 
proach could provide novel insights to the nonperturbative connection between the 
center and chiral properties. 
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I. INTRODUCTION 

Strong dynamics based on Quantum ChromoDynamics (QCD) induces two main nonper- 
turbative phenomena, “quark confinement” and “spontaneous chiral symmetry breaking”, 
which are the main themes of research on strong-coupling dynamics. While both of the 
two phenomena are caused by infrared physics based on the asymptotic freedom, it is not 
fully understood how they are related with each other. Some of lattice QCD simulations 
show that these two phenomena are almost simultaneously turned off at a certain tempera¬ 
ture [1-3], where the transition to the quark-gluon plasma (QGP) phase occurs. Regarding 
the connection between the two phenomena, we here have several questions: Do the tran¬ 
sitions take place at the same temperature? Does one of the transitions induce the other? 
Investigating these questions could lead to understanding of the nonperturbative relation 
between confining and chiral symmetry breaking. There are several approaches for study on 
the linkage between these phenomena, including the analysis on quark confinement in terms 
of chiral symmetry breaking [4-11] or vise versa [12, 13] from lattice QCD, and the effective 
models motivated by both chiral and center symmetries [14-18]. 

In the presence of exact Z 3 center symmetry, one can study the confining/deconfining 
phase transition by use of the expectation value of the Polyakov loop. Since it is the order 
parameter of center symmetry, it works as an indicator of confining/deconfining through 
the quark free energy. Quenched QCD (without dynamical quarks) [1] and adjoint QCD 
(with adjoint quarks instead of fundamental ones) [19] are examples of this case. However, 
their results on coincidence of the transition temperatures are inconsistent: For the quenched 
QCD it is reported that the transition temperatures coincides while they do not for the 
adjoint QCD. In the first place, these two theories are too far from the realistic QCD in a 
sense that one contains no dynamical quarks while the other contains an excessive amount 
of dynamical degrees of freedom of quarks. 

For the physical Nf = 2 -|- 1 QCD, where u- and d-quark masses are degenerated while 
s-quark has a heavier mass, things get more involved: the dynamical fundamental quarks 
explicitly break center symmetry, so what we can do the best is study the crossover tran- 
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sition based on the expectation values of the Polyakov loop. In the literature [2, 3], one 
hnds inconsistent results on coincidence of the crossover temperatures of center and chiral 
transitions, depending on the scheme and setup of the lattice simulations: Some works argue 
that the two critical temperatures are very close [ 2 ] while the others show that they are not 
[3], 

One ideal situation for studying this topic in details is that the exact center symme¬ 
try and the dynamical fundamental quarks are reconciled. As well-known, the existence of 
dynamical fundamental quarks breaks center symmetry explicitly, thus the above accom¬ 
modation seems infeasible. However, appropriate boundary conditions for quarks enable us 
to realize such a situation. By imposing three different twisted boundary conditions on the 
three fundamental quarks (shifted by 27r/3) in the compact imaginary-time direction, we 
realize center symmetric SU{3) gauge theory with three dynamical quarks in fundamental 
representation on x (Z 3 -QCD model) [20-24]. By investigating this model, we could 
make progress in elucidating the connection between conhning/deconhning and chiral tran¬ 
sitions. The analytical study on this theory based on the Polyakov-loop extended chiral 
model [15, 16, 18] was initiated in Ref. [23], which shows that the center symmetry is spon¬ 
taneously broken at certain temperature, associated by the manifestation of flavor symmetry 
breaking. Although the chiral model helps understand rough picture of the nonperturbative 
properties, we cannot fully eliminate the model artifacts. 

In this paper, we numerically investigate hnite-temperature Z 3 -QCD model on the lat¬ 
tice, with emphasis on center phase transition and its influence on the chiral properties. 
We formulate center symmetric SU (3) lattice gauge theory with three fundamental Wilson 
quarks by twisting the boundary conditions in the imaginary time direction. We calculate 
the expectation value of Polyakov loop and the chiral condensate on 16^ x 4 and 20^ x 4 
lattices with mps/my = 0.70 fixed. We hnd out the first-order center phase transition at 
certain temperature, where the hysteresis exists depending on the initial condition for the 
conhguration generation. We show that the chiral condensate rapidly decreases around the 
center critical temperature and has the hysteresis too. We also verify manifestation of the 
flavor symmetry breaking in the flavor-diagonal meson sector in the high-temperature phase, 
which was predicted in the effective chiral model study. 

The structure of the paper follows: In Sec. II, we briefly review the construction of Z 3 - 
QCD model in the continuum theory and introduce the expected properties based on the 
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chiral model study. In Sec. Ill, we formulate Z 3 -QCD model on the lattice and note details 
of the simulation setup including the parameter fixing for finite-temperature simulations. In 
Sec. IV, we show our results for the phase transitions associated with the center and chiral 
symmetries. Sec. V is devoted to summary and discussion. 


II. Z3-QCD MODEL 


A. Formulation as a continuum theory 


We give a brief review of Z 3 center-symmetric SU (3) gauge theory with three fundamental 
quarks on x proposed in Refs. [20]. Here, the compact dimension can be seen as 
imaginary time direction in the present work. We first consider the partition function Z in 
Enclidian spacetime with one compact dimension as 


Z= D^D^DAexp[-S], 


S = Rx 


dr 


. / ^ . 


( 1 ) 

( 2 ) 


where 4// is a quark field for degenerate three-flavor quarks with flavor index / = 1,2,3. 
Du = dfj_ — iA^ is a covariant derivative with the SU (3) gange field A^ and the field strength 
is given by = d^A^ - d^A^ - A^]. 

We here regard a compact dimension (r) as imaginary time direction, and impose an 
anti-periodic boundary condition on a qnark field in the direction as 


T = 13) = r = 0). (3) 

By adopting this bonndary condition, we now work on the finite-temperatnre system with 

/S = 1/T. 

It is known that the above action itself is invariant under the Z^ center transformation, 
which is generated by the center elements of global color SU{3) transformation. However, 
this transformation eventnally gives rise to 27r/3 shift of phase of the quark boundary con¬ 
dition as 

t = /3) = r = 0), (4) 

with k = 0,1, 2. Thereby, the Z^ center symmetry is explicitly broken via the fundamental- 
qnark bonndary conditions (3) (4) in the usnal three-flavor SU (3) gauge theory, and of course 
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in the realistic Nj = 2 + 1 QCD. 

We now consider a case that the three flavors have distinct boundary conditions as fol¬ 
lowing, 

^f{x, t = /3) = r = 0), (5) 

for / = 1,2,3. It is nothing but SU{3) gauge theory in the presence of three fundamental 
quarks with flavor-dependent boundary conditions. 

For this case, the transformation in Eq. (4) shifts the quark boundary conditions as 

^/(F, t = /3) = r = 0). (6) 

The twisted angles for the three flavors are changed into k + f — 1, but this can be straight¬ 
forwardly relabeled as / — 1 and returns back to the first place in Eq. (5). In other words, 
we can rename the flavor f + k as f (mod 3). This means that the three-flavor SU{3) gauge 
theory with the special twisted boundary condition Eq. (5) is invariant under the center 
transformation. We call this exactly-center-symmetric model as Z^-QCD model [20-24]. We 
note that the flavor-dependent twisted boundary condition is translated into the insertion of 
the flavor-dependent imaginary chemical potential by use of gauge transformation as shown 
in Appendix. A. Hereafter, we use f = u,d,s instead of / = 1, 2, 3 as indices for flavor. 

B. Chiral structure and flavor symmetry in Z3-QCD model 

We here comment on flavor symmetry possessed by Z3-QCD model in the massless limit. 
In Z3-QCD, due to the boundary conditions, the flavor-chiral SU{3)l x SU{3)ji symmetry 
is explicitly broken to its Cartan subgroup x generated by A3, Ag of the Gell- 

mann matrix elements. However, this influence from the boundary condition disappears in 
the infinite compact-circumference limit, or in the zero-temperature limit. Thus, in this 
limit, Z3-QCD is reduced to the standard three-flavor QCD at zero temperature. 

Eor nonzero temperature, the effects of the twisted boundary condition exists in principle, 
namely the flavor-chiral symmetry is broken to C(l)| x C(l)|. at the action level. In this 
work, we regard this Cartan subgroup as the specific flavor-chiral symmetry of the Z 3 -QCD 
model. Thus, the pattern of the chiral symmetry breaking is basically given by C(l)| x 
uii)l —)■ in this model. While the chiral transition in the standard three-flavor QCD 
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model is believed to be first-order in the chiral limit based on the universality class [25], the 
same discussion is probably not valid for the Z 3 -QCD model. 

Here, we make a supplemental comment on the flavor-chiral symmetry 17(1)1 x U{1)\ 
in the Z3-QCD model discussed above. It is notable that, indeed, we do not know at 
how high temperature the physical quantities start to be affected by the twisted boundary 
condition. It may be as soon as the temperature is turned on or may be at a certain nonzero 
temperature. Actually, as we will see in the next subsection, the study based on the chiral 
effective model indicates that the chiral condensate is insensitive to the boundary condition 
below the center phase transition temperature, where the flavor symmetry breaking due to 
the boundary condition is not manifest in the effective thermodynamic potential. 

The signihcance of this model is that we can study the phase transition with respect 
to the Z 3 center symmetry by calculating expectation values of Polyakov loop, even in 
the presence of the dynamical fundamental quarks. We here have physical and theoretical 
questions: If the phase transition occurs, what is the order of the phase transition? How is 
the chiral condensate affected by the center phase transition? We study these topics in the 
next section numerically, although our simulations are not carried out in the chiral limit. 
The study on such an ideal model with the exact center symmetry may help elucidate the 
relation of center and chiral properties from novel viewpoints. 

C. Expected properties from chiral model and QCD with a finite chemical potential 

In the study of Z3-QCD based on the Polyakov-loop-extended NJL (PNJL) model [15, 16, 
18], the following three properties are predicted [23] (we call this model just as the effective 
chiral model in the present paper.): 

(i) Z 3 center symmetry is spontaneously broken in the high-temperature phase, where 
the order of the phase transition is first. 

(ii) Although the SU (3) flavor symmetry the standard three-flavor QCD has is spoiled by 
the twisted boundary condition at the action level, the effective thermodynamical potential 
of the PNJL model in the low-temperature phase is not affected by the boundary condition. 
It is expected that the Z3-QCD model becomes SU{3) flavor-symmetric in the phase. On 
the other hand, the flavor symmetry breaking becomes manifest in the high-temperature 
phase. 
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(iii) In a chiral limit, even above the critical temperature of the center symmetric phase 
transition, the chiral condensate has a nonzero value and the chiral symmetry is still broken. 
It might be an artifact coming from the model cutoff in the chiral effective model. Instead, 
the value of chiral condensate has a specific jump at the center critical temperature. 

The present model is also related to the works on QCD with hnite imaginary chemical 
potential. As discussed in Refs. [26, 27], the partition function has 27r/3 periodicity in 
the imaginary chemical potential. Thus we speculate that the critical temperature for the 
chiral phase transition in the present model, which also has 27r/3 periodicity of the temporal 
direction, could be the same as that of the standard three-flavor QCD. 


III. SIMULATION SETUP 


We utilize the Iwasaki gauge action with naive Wilson fermions in our lattice numerical 
simulation. The dehnition of the action is given by 


S 

S, 


s, + Sf, 


( 7 ) 


( 8 ) 

f=u,d,s x,y 

where {3 = Q/g"^, in which g is a lattice bare gauge coupling constant, ci = —0.331, cq = 1— 8 ci 
and and denote the plaquette and rectangular, respectively. In Eq. ( 8 ), 

4 

^x,y — 3x,y — ^ ^ ^ {(7 ~ 'ly)Ux,y3x+jl,y + (1 + 1 y)Uy ^5x^y+■ ( 9 ) 

y=l 


Here k is the hopping parameter. We note that the value of k in Z 3 -QCD model is universal 
for all flavors. 

To realize the twisted boundary condition in Eq. (5) on the lattice, we introduce the 
following boundary conditions for the link variable only in the fermion action Eq. ( 8 ): 


Ui{x,T = Nr) = —Ui{x,T = 0 ) for M-flavor , 
774 ( 0 ;, r = Nr) = 7 / 4 ( 0 ;, r = 0) for d-flavor, 

Ui{x,T = Nr) = — 7 / 4 ( 0 ;, r = 0 ) for s-flavor. 


( 10 ) 



These conditions from top to bottom are the same with those for the standard hnite- 
temperature QCD with imaginary chemical potential; fij = 0,27r/3 and 47r/3, respectively. 
We have to use the Rational Hybrid Monte Carlo (RHMC) algorithm to calculate the fermion 
determinant for each flavor, since three fermions have a different boundary condition with 
each other. 

Firstly, we perform the zero-temperature simulation using 16^ lattices to obtain the line 
of constant-physics. We carry out the simulations for several values of hopping parameter 
K, with each {3 value and measure the flavor-singlet pseudo-scalar mass {mps) and vector 
meson mass {my) for each flavor. The number of trajectories we generate is 2, 000 - 3, 000, 
and we measure the correlator of these hadronic states every 10 Monte Carlo trajectories. 
The estimated autocorrelation length is around 100 trajectories in this simulation. We 
summarize results on mass measurement in Table II and III in Appendix C. 
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1.30 

1.40 

1.50 

1.55 

1.60 

1.70 

1.80 

1.90 

1.95 

2.00 

2.10 

2.20 

K 

0.2019 

0.1975 

0.1921 

0.1892 

0.1861 

0.1793 

0.1725 

0.1663 

0.1636 

0.1611 

0.1571 

0.1539 


TABLE I: Simulation parameters: /I ^.nd k 


We £x the ratio between mps and my constant, namely mps/my = 0.70, and tune 
the value of hopping parameter for each (3 shown in Table I. Using these parameter sets in 
Table I, we perform the hnite temperature simulation on 16^ x 4 and 20^ x 4. The number of 
trajectories for the hnite temperature lattice setup is 500 - 5, 000. We measure the Polyakov 
loop in temporal direction for every Monte Carlo trajectory and the chiral condensate every 
10 trajectories. 

We note that the masses of havor-singlet mesons are the same with that of the standard 
three-havor QCD in zero-temperature. It shows that, at least in the havor-singlet sector, 
the breaking of SU{3)l x SU{3)p havor-chiral symmetry due to the Z 3 twisted boundary 
condition (5) is not observed for zero-temperature. It is consistent with the general argument 
in the previous section and the result of PNJL model in Sec. IIC. 
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FIG. 1: Polyakov loop distribution plot in Z 3 -QCD (left) and the standard three-flavor QCD 
(right). Based on 16^ x 4 lattice for /3 = 1.70, 2.00, 2.20 with the same values of k in both panels. 

IV. SIMULATION RESULTS 


A. Polyakov loops and center symmetry 

We first show the existence of center symmetry in the present model based on the 
distribution plot of Polyakov loop (L), 

'N^ 


r 1 V- 1 

L = — y - tr 
U 4^ 3 

X 


Y[Ur{x,i) 


2 = 1 


( 11 ) 


Here V denotes the spacial volume in a lattice unit. As shown in the left panel of Fig. 1, 
the Polyakov loops are distributed around the origin in the low (3 regime while three vacua 
exist in the high jS regime for Z3-QCD. On the other hand, those in the standard three- 
flavor SU{3) gauge theory in the right panel of Fig. 1 indicate explicit breaking of Z 3 
center symmetry. These results obviously show that the Z3-QCD model possesses exact Z 3 
center symmetry at the action level while it seems to undergo spontaneous breaking of the 
symmetry in the high-temperature phase. 

Next, we investigate temperature dependence of the Polyakov loop by varying (3 along 
with the line of constant-physics, namely mpg/my = 0.70, shown in Table I. We generate 
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FIG. 2: 13 dependence of the magnitude of Polyakov loop ((|F|)) for the Z 3 -QCD and standard 

three-flavor QCD on 16^ x 4 (left) and 20^ x 4 (right) lattices. For the Z 3 -QCD model, the data of 
(|L|) started with the cold start (triangle (blue) symbols) have a clear jump from zero to non-zero 
values around the region 1.55 < /3 < 1.60 in both panels, while the jump occurs in 1.80 < /3 < 1.90 
(left) and 1.90 < (3 < 1.95 (right) for the data generated by the hot start (circle (red) symbols). In 
the regions between these two jumps, the hysteresis exists in Z 3 -QCD model. On the other hand, 
the data of the standard three-flavor {Nj = 3) QCD (square (black) symbols) do not show such a 
jump from zero to nonzero nor hysteresis. 


configurations with two types of initial condition; cold start and hot start. In both panels 
of Fig. 2, the triangle (blue) symbol denotes the data started with cold start. The corre¬ 
sponding initial conhguration lives in the ordered phase, and we set all initial link variables 
to unity. On the other hand, the circle (red) symbol denotes the ones started with hot start. 
The corresponding configuration is in the disordered phase, and the initial link variable 
is a random number. The square (black) symbol shows the result of the standard three- 
flavor QCD with the periodic boundary condition for spacial directions and the anti-periodic 
condition for temporal direction with the same values of [3 and k as Z 3 -QCD simulations. 

Now, let us look into the results in details. 

Firstly, for the Z 3 -QCD model, we find hysteresis in the range of 1.55 < (3 < 1.90 
depending on the initial conditions (cold or hot). On the other hand, we hnd that there are 
no hysteresis in the data for the standard three-flavor QCD. We note that the hysteresis is 
a signal of the hrst order phase transition. 

Secondly, in the low-temperature phase, the magnitude of Polyakov loop is exactly zero for 
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p 


FIG. 3: Polyakov loop susceptibility X(|L|> as a function of /3 on 20^ x 4 lattices. Circle (red) 
and triangle (blue) symbols denote the data generated with hot and cold starts in Z 3 -QCD model, 
respectively. Square (black) symbol shows the data for standard three-flavor [Nf = 3) QCD model. 

the Z3-QCD model. It originates in the existence of exact Z3 center symmetry. From these 
results, we argue that the Z 3 -QCD model undergoes hrst-order phase transition, where the 
Z 3 center symmetry is spontaneously broken, while the standard three-flavor QCD undergoes 
the crossover transition. 

Figure 3 shows the Polyakov loop susceptibility dehned by 

X(|i|) = r [(|L|^> - {|L|)"] . (12) 

The meaning of colors of symbols is the same as the one in Fig. 2. The data for the Z 3 - 
QCD model has a relatively clear signal of peaks around [3 = 1.60 and [3 = 1.90 for cold 
and hot starts, respectively. The peak appears because of the co-existing states between 
the two phases associated with the phase transition [28]. On the other hand, the standard 
three-flavor QCD does not show a clear transition point. Actually, such a situation makes 
it difficult to determine the critical temperature in the physical QCD [29-31]. 

We also investigate the volume dependence of the critical (3 and hnd small hnite volume 
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effect, which we will not discuss further in the present paper. The precise determination of 
the critical temperature in the large volume and continuum limits remains as a future work. 

Our results on center symmetric phase transition indicates the following points: 

(1) Z3-QCD model on the lattice possesses exact Z 3 center symmetry. 

( 2 ) Z 3 center symmetry is spontaneously broken in the high-temperature phase. 

(3) The order of the center phase transition is hrst. 

B. Chiral condensates 

Here, we investigate the chiral property in the hnite-temperature Z3-QCD mode, which is 
characterized by the chiral condensate. In the chiral limit, the non-vanishing chiral conden¬ 
sate is a signal of the spontaneous chiral symmetry breaking. As we discussed in Sec. IIB, 
because of the twisted boundary condition, the symmetry breaking pattern is expected to 
be Uil)l X Uil)l —)■ in the present model. On the other hand, as we have discussed 

in Sec. IIC, the result of the chiral model implies that the flavor symmetry breaking due to 
the flavor-dependent twisted boundary condition gets manifest only in the high-temperature 
phase [23]. To investigate this conjecture, we will investigate the chiral condensate for each 
flavor independently. 

Since our simulation is performed for heavy mass region mps/my = 0.70 and utilize 
the Wilson fermion, no chiral phase transition takes place. Nevertheless, we expect that, 
in the system with massive fermions, decrease of chiral condensate indicates effective and 
approximate restoration of chiral symmetry. Actually, several lattice numerical simulations 
in these years reveal the chiral property for the realistic Nf = 2-1-1 QCD at hnite temperature 
based on chiral condensate and its susceptibility as the effective order parameter. We note 
that the recent lattice results using the staggered fermions and domain-wall fermions give 
the consistent chiral critical temperature [29-31] . 

We consider the following flavor-diagonal expectation value of the subtracted chiral con¬ 
densate [32-34] for each flavor [35-37] (we do not consider flavor-mixing condensate in the 
present paper.), 

= (2mpcAc)(2K)^^(P(T,t)P'^(0,0)). 


X 


(13) 
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Here, rupcAc is partially conserved axial current (PCAC) mass, and P denotes the pseudo 
scalar state defined by P = for the corresponding flavor (/) in LHS. The PCAC 

mass is defined via axial Ward identity; 


2mpcAc 


E.-^4(A4(T,t)Pt(0,0)) 

E.^(P(pt)pt(0,0)) 


(14) 


where corresponds to the axial vector current defined by Here again, 

the label of flavor (/) is fixed. The values of mpcAc are summarized in Appendix C. 

Before showing our results of numerical simulation for chiral condensates for each flavor, 
we comment on the notation of flavor in this study. To fix the name of quarks, we firstly 
observe the complex phase (0) of Polyakov loop given as L = |L|e*'^ for each configuration. 
As discussed in the previous section, in the center-symmetric phase the definition of the 
complex phase is meaningless since (|T|) = 0 , while in the center broken phase we can define 
the value of 0 [46]. Next, we measure the correlators in Eq. (13) using the Dirac operator 
given in Eq. (9), where the link valuable has the boundary condition given by 


Ur{x, T = N^) = -P^Urix, 0), (15) 

where 6 takes value 0 or ±27r/3. We define the flavor u as in 

(j) + 9 = 0 (mod 27r) for w-fiavor, (16) 

while the flavor d and s are defined with the following total phases; 

0-1-6' = 2tt/3 (mod 27r) for d-fiavor, (17) 

(j) + 9 = Att/S (mod 27r) for s-fiavor. (18) 


Figure 4 shows the /3 dependence of the chiral condensate for m- flavor of the Z^- and three- 
flavor QCD. Circle (red), triangle (blue) and square (black) symbols stand for those data 
associated with the hot start and cold start in Z 3 -QCD model and the standard three-flavor 
QCD at finite temperature, respectively. All results have common qualitative properties, 
where the chiral condensate gradually decreases as 0 increases. It is also notable that P 
dependence of chiral condensates in Z 3 -QCD model has hysteresis between the cold and hot 
starts as with that of the Polyakov loop. These results imply that effective restoration of 
the chiral symmetry is progressing from 0 = 1.55 to 0 = 1.95 for the Z 3 -QCD model. These 
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FIG. 4: /3 dependence of the expectation values of subtracted chiral condensates for 

■^ 3 -QCD and three-flavor QCD on 20^ x 4 lattices. Circle (red), triangle (blue) and square (black) 
symbols stand for those data associated with the hot start and cold start in Z 3 -QCD model and 
the standard three-flavor QCD, respectively. The error bar denotes the statistical error coming 
from the pseudo-scalar correlator in Eq.(13). 

values of (3 where the hysteresis exists are the same as those of the center phase transition. 
Furthermore, the decreasing rate of chiral condensate in Z3-QCD model is relatively larger 
than the one in the standard three-flavor QCD. 

According to the arguments in Refs. [26, 27], the partition function for QCD has 27r/3 
periodicity in the imaginary chemical potential, and it is expected that the chiral phase 
transition temperature in our model is the same as that of the standard three-flavor QCD. 
On the other hand, since the flavor-chiral symmetry in Z 3 -QCD model in the chiral limit is 
broken to U{1)\ x D(l)|j. due to the twisted boundary condition, we may have the smaller 
number of Nambu-Goldstone modes than the usual three-flavor QCD, which lifts the phase 
transition temperature in general [38]. (As we have discussed, the chiral effective model 
indicates that the full flavor symmetry is effectively preserved in the low-temperature phase, 
and in such a case we should have the common number of Nambu-Goldstone modes.) Our 
results in Fig. 4 indicate qualitatively the same chiral crossover temperature in Z 3 -QGD 
and three-flavor QGD, while, to be quantitative, the temperature in Z 3 -QGD seems slightly 
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/3 


FIG. 5: /3 and flavor dependences of the expectation values of subtracted chiral condensates (qq) 
for each flavor in Z 3 -QCD model. The lattice size is 20^ x 4. Circle (red), square (orange) and 
triangle (violet) symbols denote u-, d- and s-flavor generated with hot start, respectively. 

higher than that in the three-flavor QCD. We do not yet have sufficient ingredients to con¬ 
clude on this question. Higher statistics and investigation of its susceptibility are necessary 
to determine the critical temperature in the massless limit. It is of well-known difficulty to 
determine the critical temperature of the chiral phase transition [39]. 

We also note that the absolute values of the chiral condensates in Z 3 - and three-flavor 
QCD are different in the low-temperature phase, which may indicate the qualitative differ¬ 
ence of the chiral property between the two theories. 

Next, we focus on the flavor symmetry breaking in the high-temperature phase. Figure 5 
shows the expectation values of chiral condensates for each flavor. Here, circle (red), square 
(orange) and triangle (violet) symbols denote m-, d- and s-flavor generated with hot start, 
respectively. Three components of chiral condensate are degenerate in the low-temperature 
phase. On the other hand, in the high-temperature phase, there appears clear flavor symme¬ 
try breaking. Two of them, whose total complex phase {4>+0) are nontrivial, are degenerated 
because of the momentum shift of the twisted boundary condition. It indicates that at least 
the Z 3 center of SU{3) flavor symmetry, which commutes with the Cartan subgroup, is ef¬ 
fectively preserved in the low-temperature phase, while the breaking of this symmetry gets 
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manifest in the high-temperature phase. Although our simulation cannot fully verify the 
conjecture of the chiral model shown in (ii) of Sec. IIC, which states that the SU (3) flavor 
symmetry is not affected by the twisted boundary condition in the low-temperature phase, 
the above result on the flavor symmetry is consistent with this conjecture. On the other 
hand, Fig. 4 shows the different values of chiral condensates between the and three-flavor 
QCD, thus the two theories seem to have difference in chiral and flavor properties. Further 
study is required to understand the structure of flavor and chiral symmetry in the Z 3 -QCD 
model. 


V. SUMMARY AND DISCUSSION 

In this work we numerically simulate the hnite-temperature Z 3 -QCD model on the lattice 
by introducing flavor-dependent twisted boundary conditions, with emphasis on center phase 
transition and its influence on the chiral properties. We calculate the hnite-temperature 
expectation value of Polyakov loop and the havor-diagonal chiral condensates on 16^ x 4 
and 20^ x 4 lattices on the line of constant physics realizing mps/my = 0.70. We hnd out 
the hrst-order center phase transition at a certain temperature with the hysteresis depending 
on the initial conditions for conhguration generation. The chiral crossover transition takes 
place around the critical temperature of center transition, where it also has hysteresis in the 
same range as that of the center phase transition. By comparing this approximate chiral 
restoration to that of the standard three-havor QCD, we hnd that the decrease of chiral 
condensate in Z 3 -QCD is more rapid while temperatures of the two crossover transitions 
are almost the same. We also obtain the result which supports the manifestation of havor 
symmetry breaking due to the boundary condition in the high-temperature phase, which 
was predicted in the study of the ehective chiral model. 

One of motivations for this work is study how the chiral condensate reacts to the center 
phase transition. Our preliminary result indicates not only that the chiral condensate has 
rapid decrease around the center critical temperature, but also that it has the hysteresis in 
the same range as that of the center phase transition. Based on this result, we argue that 
the chiral and center properties at hnite temperature have a strong correlation at least in the 
present model. One possible reason for the correlation is that the hrst-order center phase 
transition might work as a trigger to the rapid decrease of chiral condensate. On the other 
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hand, the argument in Ref. [27] implies that the chiral critical temperature in the present 
model is the same as that of the standard three-flavor QCD, thanks to the 27r/3 periodicity 
of imaginary chemical potential. Our result on comparison of the two chiral (crossover) 
transition temperatures is consistent with this conjecture. The precise study on the critical 
temperature in the Z 3 -QCD model may help determine chiral crossover temperature Tc in 
the standard three-flavor QCD. 

By comparing our results to those of the effective chiral model, we have supported the 
several conjectures based on the chiral model; first-order center phase transition, latency of 
flavor symmetry breaking in the low-temperature phase, and manifestation of flavor symmetry 
breaking in the high-temperature phase. In particular, we have observed that Z 3 center part 
of the SU{3) flavor symmetry seems to be intact in the low-temperature phase, while its 
breaking becomes manifest in the high temperature phase, accompanying the spontaneous 
breaking of center of the SU{3) color symmetry. This is consistent with the conjecture 
of the effective chiral model, but is not sufficient to draw a conclusion. On the other hand, 
we could not hnd the specihc jump of the chiral condensate right at the center critical 
temperature, which was seen in the chiral model as discussed in (iii) of Sec. IIC. Although 
our simulation is not performed in the chiral limit, we may be able to interpret the jump 
seen in the chiral model as a model artifact and a remnant of the rapid decreasing of the 
chiral condensate seen in our simulation since the cut-off effect may be visible at the high- 
temperature in the effective model. Further study is required to clarify the flavor and chiral 
properties in the Z 3 -QCD model. 


For future works, we address following points. 

Towards a chiral limit: By approaching a chiral limit, the Z 3 -QCD model has both exact 
center and exact chiral symmetries in presence of fundamental quarks. We may be able to 
determine the both center and chiral critical temperatures using the exact order parameters, 
Polyakov loop and chiral condensate. Therein, we can compare the two critical temperatures 
and discuss relation of conhning and chiral properties in more details. 

Towards smaller lattice and larger volumes: To obtain the critical temperatures, we have 
to take a continuum limit to remove a lattice artifact and take large volume limit to estimate 
hnite volume effects. Several works for the three-flavor real QCD near the physical points 
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have been progressing, so that in principle it is doable at present. Comparison with the 
critical temperatures between the Z 3 -QCD and the standard three-flavor QCD in these 
limits must be interesting for understanding the center and chiral structure. 

Towards topological objects (fractional instantons): In the presence of center-symmetric 
Polyakov loop holonomy, or equivalently the center-symmetric twisted boundary conditions, 
fundamental topological objects become fractional instantons with fractional topological 
charge as 1/N with N being the rank of gauge group [40, 41]. In our simulations, we expect 
that 1/3 fractional instantons are present, and have influence on physical observables such 
as topological susceptibility. It could be valuable to investigate the topological properties in 
relation to the recent interest in the novel topological objects [42-44]. 

Towards further application of gauge theory with twisted boundary conditions: Our result 
is of significance in a sense that we first observed the center first-order phase transition in 
the lattice QCD simulation with dynamical fundamental quarks by introducing the twisted 
boundary condition (See also Appendix B). It is notable that the only difference between 
the standard three-flavor QCD and our Z3-QCD model is the twisted boundary condition^ 
and this difference disappears in the zero-temperature limit. It means that, by choosing 
appropriate boundary conditions for quarks, we can realize an ideal situation that the con- 
flning/deconflning properties are well described by center symmetry and Polyakov loop. 
We consider that such a technique using the twisted boundary conditions helps understand 
broader topics on nonperturbative aspects of strong-coupling physics. 
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Appendix A: The flavor-dependent b.c. as imaginary chemical potential 

Consider the transformation where the fermion helds in Eq. (5) are transformed as 

. (Al) 

Then, the action (2) is translated into 

S = J d^x dr ^'E/( 7 ^T)^ + i(/-l) 74 // 3 -hm)T/ + . (A 2 ) 

This is the theory with the flavor-dependent imaginary chemical potentials and the usual 
anti-periodic boundary conditions Eq. (3). 

Appendix B: Comment on the other possible Z3 symmetric theory 

In the canonical ensemble, it is known that the Polyakov loop has exact zero expectation 
value [45], which is similar to that of Z 3 -QCD. Here, we briefly discuss the difference between 
the canonical ensemble and our Z 3 -QCD model. 

The grand canonical partition function can be written by 

Zgc{T, fi) = J det M{U, p), (Bl) 

where jj, denotes a quark (real) chemical potential. Using ^GC the canonical partition 

function can be expressed as 

/ OO 

d (^) e-^i^'/^ZacX, ii = t^„). (B2) 

with the quark number Q. 

By using (27rT/3)-periodicity of the grand canonical partition function as ZQc(T,i{^i + 
2nT/3)) = ZGciT,i^j), the canonical partition function can be expressed as 

" ^ (t ) V/) 


Zc{T,B) 


(B3) 
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for the baryon number B{= Q/3), apart from a normalization factor. 

Insertion of the imaginary chemical potential with 

fii = ^— for A; = 0,1,2, (B4) 

can be translated by the center transformation {z{k) = acting only on the link 

variables in the fermion action, 

U4{x,X4 = 0) ^ z{k)U4{x,X4 = 0). (B5) 

Therefore, the Dirac determinant satishes 

det M{z{k)U, ifii) = det M(17, ip/ — i27rTk/3). (B 6 ) 

Using the periodicity of ZQciT^ijjij) and the above relation, the canonical partition func¬ 
tion is expressed as 

Zc(T,B) = i (^) 

+ ZQQ{T,i^! — z2;rT/3) + ZQQ[T^iji] — z4;rT/3)] , 

2 

^ dgt M{z{k)U, ifii), (B7) 

k=0 

which means that Zc{T,B) is the average of three center sector. It leads the exact zero 
expectation values of the Polyakov loop as 

oc 1 + = 0. (B 8 ) 

On the other hand, the Z 3 -QCD model introduces the flavor dependent imaginary chem¬ 
ical potential. Thus, the partition function is given by 

P 2 

^Z 3 -qcd(T) = / Y[ det M(U, p = i27TkT/3). (B9) 

k=0 

This partition function keeps center symmetry, since the integrand is invariant under the 
center transformation as follows: 

det M(f/, ju = 0) det M(U, ju = i27rT/3) det M(U, p = —i27rT/3), 

^ det = 0) det p = i27rT/3) det p = -i27rT/3), 

= det M{U, jj, = —i27iT/3) det M{U, /i = 0) det M{U, p = i27iT/3). 



(BIO) 



21 




FIG. 6: The hopping parameter k dependence PS and PCAC masses. Lines denote fit results 
using (1 /k — l/Kc)-type function. 

Appendix C: PCAC mass and PS meson mass 

Here, we summarize pseudo-scalar (mps), vector (my), their ratio (mps/my), and PCAC 
(mpcAc) masses at the zero-temperature simulation using 16^ lattice in Table II and III. 
Figure 6 shows the hopping parameter dependence of mpcAC and mpg with the £t results 
using (1 /k — I/kc). For smaller masses, both are proportional to (1 /k — 1/kc), and they 
become zero at almost the same k up to 0{ampcAc), which are expected behaviors of Wilson 
type fermion [37]. For the tuning of the mass parameter at mpg/my = 0.70, we interpolate 
mps and my as a function of {1 /k — 1/kc) and determine the line of constant physics. 
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/? K mpsa my a mps/my mpcAC® 

1.80 0.1600 1.3487(36) 1.4365(50) 0.9389(41) 0.1715(3) 

1.80 0.1700 0.8583(25) 1.0399(28) 0.8254(33) 0.0888(1) 

1.80 0.1720 0.6644(26) 0.8944(37) 0.7428(42) 0.0565(1) 

1.80 0.1725 0.5906(26) 0.8392(40) 0.7038(46) 0.0459(1) 

1.90 0.1600 1.0644(49) 1.1807(60) 0.9015(62) 0.1233(6) 

1.90 0.1650 0.6968(34) 0.8849(51) 0.7874(59) 0.0663(1) 

1.90 0.1660 0.5710(43) 0.7905(43) 0.7223(67) 0.0510(1) 

1.90 0.1665 0.4920(31) 0.7310(52) 0.6731(64) 0.0367(1) 

1.95 0.1610 0.7852(28) 0.9393(35) 0.8359(43) 0.0866(1) 

1.95 0.1620 0.6959(29) 0.8642(39) 0.8053(49) 0.0731(1) 

1.95 0.1630 0.5861(31) 0.7871(40) 0.7446(55) 0.0526(1) 

1.95 0.1640 0.4325(47) 0.6754(51) 0.6404(85) 0.0351(1) 

2.00 0.1550 0.9948(48) 1.0942(91) 0.9092(87) 0.1167(23) 

2.00 0.1600 0.6203(28) 0.7914(42) 0.7838(55) 0.0633(1) 

2.00 0.1610 0.4948(31) 0.6996(35) 0.7073(57) 0.0457(1) 

2.00 0.1615 0.4231(37) 0.6524(44) 0.6485(72) 0.0299(1) 

2.10 0.1525 0.8455(25) 0.9464(36) 0.8934(43) 0.1180(7) 

2.10 0.1550 0.6563(16) 0.7934(27) 0.8272(35) 0.0774(1) 

2.10 0.1560 0.5683(24) 0.7244(33) 0.7845(49) 0.0623(1) 

2.10 0.1570 0.4576(25) 0.6423(34) 0.7124(54) 0.0406(1) 

2.20 0.1500 0.7675(31) 0.8677(31) 0.8845(48) 0.1040(3) 

2.20 0.1525 0.5769(27) 0.7107(30) 0.8117(51) 0.0636(1) 

2.20 0.1540 0.4354(25) 0.6283(41) 0.6930(60) 0.0383(0) 

2.20 0.1545 0.3845(41) 0.6116(50) 0.6287(84) 0.0285(1) 

TABLE III: Summary of pseudo-scalar, vector, and PCAC masses for (3 = 1.80 — 2.20 using 16^ 

lattice. 



